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Ratios of quadratic forms in correlated normal variables which introduce noncentrality into the 
quadratic forms are considered. The denominator is assumed to be positive (with probability 
1). Various serial correlation estimates such as least-squares, Yule- Walker and Burg, as well 
as Durbin- Watson statistics, provide important examples of such ratios. The cumulative dis- 
tribution function (c.d.f.) and density for such ratios admit saddlepoint approximations. These 
approximations are shown to preserve uniformity of relative error over the entire range of sup- 
port. Furthermore, explicit values for the limiting relative errors at the extreme edges of support 
are derived. 
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1. Introduction 

Consider the ratio of quadratic forms 

where, without loss in generality, A and B are assumed to be n x n symmetric. Let 
N(/x, I n ) and suppose that B is also positive semidcfinite, thereby ensuring that the 
denominator is positive with probability one. There is no loss in generality in having the 
covariance of e as the identity. This is because, if the distribution of e were N(/x, £), 
then (1) describes the model with S 1/2 AS 1/2 and S 1/2 BS 1/2 replacing A and B, 

1/2 

respectively, and £ ' fi replacing /j, in the distributional assumption on e. Thus, model 
(1) incorporates all dependence among the components of e, as well as noncentrality that 
occurs when 0. 

Various types of saddlepoint approximations for the distribution (c.d.f.) and density 
of R have been proposed, beginning with the seminal work on serial correlations in 
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Daniels (1956). Further marginal distributional approximations are given in McGregor 
(1960), Phillips (1978), Jensen (1988), Wang (1992), Lieberman (1994a, 1994b) and 
Marsh (1998). Joint distributional approximations for the set of serial correlations com- 
prising the correlogram were initiated by Daniels (1956) and continued in Durbin (1980) 
and Butler and Paolella (1998). 

The main contributions of the current paper are in establishing the uniformity of 
relative errors for the saddlepoint cdf and density approximations in the right tail when 
used with univariate ratios R in a class C-r, that is defined below. This class encompasses 
all of the examples in the aforementioned papers. Expressions for the limiting relative 
error are given as the right edge of support for R is approached and the sample size n 
is held fixed. These expressions are explicit in the more elementary settings in which a 
certain defining eigenvalue is simple and mostly implicitly defined when it is multiple. 
The noncentral beta distribution provides an important example in which the limiting 
error is explicit, but the defining eigenvalue is multiple. 

The left tail of R is dealt with by changing A to —A, thus switching the left tail of 
R to the right tail of —R. The results for the right tail of — R can now provide similar 
uniformity results for the left tail, when applicable. If —R is a member of class C-jz, then 
we say that R is in Cc- However, for the most part, this paper concentrates on the class 
C n - 

The class of ratios Cn is technically characterized in terms of a sequence of largest 
eigenvalues. Let (I, r) be the support of R with r possibly infinite and define X n (r) as the 
largest eigenvalue of A — rB for r <E (I, r). The class C-n is characterized as those ratios R 
whose matrices A and B admit the limit 

= limA n (r) (2) 

r — >r 

with dimension n fixed. The class C-n contains the subset B that consists of all ratios with 
bounded support (I, r); this is a property of R that is guaranteed when B > 0, or positive 
definite. Such ratios include the Durbin-Watson statistics, as well as the Yule-Walker 
and Burg estimators of serial correlation with arbitrary lag computed from least-squares 
residuals. If B > has at least one zero eigenvalue, then r may be finite or infinite. The 
portion of C-n with r = oo includes least-squares estimators of serial correlation in various 
types of models with arbitrary lag and computed from residuals with trend or covariates 
removed. Such models include those with autorcgrcssive lag in the dependent variable 
and those with lag in the additive noise. 

Large sample space asymptotics of the type considered here have not been previously 
considered for the class B. The only previous consideration for a member of the class 
C-jz — B is in Jensen (1988), (1995), Chapter 9.4. Results obtained there are in agreement 
with those below for the least-squares estimator of lag-one serial correlation when the 
time series is a mean-zero AR(1) model. 

The class Cn excludes _F-statistic and Satterthwaite-type ratios which have been con- 
sidered in Butler and Paolella (2002). In this work, A,i(?') > does not depend on r. For 
this setting, saddlepoint uniformity is also maintained; however, a different asymptotic 
saddlepoint behavior results from these different assumptions. 
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Some alternative large sample size asymptotics for the lag-one least-squares estimator, 
showing that the error is 0{n~ l ) and 0(n~ 3 / 2 ) on compact sets as n — » oo, are given 
in Lieberman (1994b) and Jensen (1995), Chapter 9.4, respectively, when ^ = 0. Such 
asymptotics, in which n — ► oo, are not considered in this paper. 

Since most of the proofs in this paper arc long and technical, they have been relegated 
to the accompanying technical report, Butler and Paolella (2007). 



2. Saddlepoint approximations 
2.1. Distribution theory 

The cdf for R in the most general setting with noncentrality is 

Pr(i? <r) = Pr ( <rj= Pr(e'(A - rB)e < 0) 



(3) 



= Pr(X r <0), 
where X r = e'(A — rB)e 

A-rB = P;A r P r , (4) 
where P r is orthogonal and A r = diag(Ai, . . . , A n ), with 

Ai=Ai(r)<---<A„ = A„(r) 

consisting of the ordered eigenvalues of (4). Whenever convenient, we suppress the de- 
pendence of the various quantities on r. The distribution of X r is therefore 

n 

X r = ^\ iX 2 (l,vt), (5) 

where {vf\ are determined as (y\, . . . , v n )' = v r = P r fi and represent the noncentrality 
parameters of the independent noncentral Xi variables specified in (5) . The ordered values 
of {Ai} are in one-to-one correspondence with the components of v r specified through 
the particular choice of P r . Notationally, we use x\ f° r the central chi-square instead of 

x 2 (M). 

Before proceeding with the development of a saddlepoint approximation for the distri- 
bution of R, we must first characterize the support of R, its relationship to the eigenvalues 
Ai(r) and A n (r) and the convergence strip for the moment generating function of X r . 

Lemma 1. All of the eigenvalues of A r are strictly decreasing in r when B > and 
decreasing when B > 0. 



Lemma 2. The distribution of R is degenerate at a single point if and only if A = cB 
for some scalar constant c. 
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A description of the support of R requires the consideration of the various cases in- 
volved which depend on eigenvalue decompositions of A and B. Suppose that B has 
p > zero eigenvalues and let O b be the orthogonal matrix of eigenvectors for B such 
that 



O b BO' b 

Denote 



A B 

Opxp 



O b AO b =(^ 11 

\<^21 ^22 

where Cn is (n — p) x [n —p) and C22 is p x p. Let iV(Ci2) denote the null space in 5i p 
for matrix C\2- 



Lemma 3. The support of R is specified in the following set of cases. 

1. Suppose that B > ; hence p = 0, and that A has rank of at least one. The support 
of R is then the finite interval (l,r) with I and r being the smallest and largest 
eigenvalues o/B _1 A. 

2- Ifp>l, so that B has at least one zero eigenvalue, then the right edge r is given 
as follows: 

(a) if C22 has a positive eigenvalue, then r = oo; 

(b) if C22 < then r < 00 and r is the largest eigenvalue o/Ag 1 (Cn — Ci2C^ 2 1 C2i); 

(c) if C22 < and C22 has at least one zero eigenvalue, then r = 00 if 
N{C22) % N(C\2); otherwise r < 00 and is the largest eigenvalue of Ag 1 (Cn — 
Ci20ciAq 1 0^, 1 C2i). Here, 0' c = (Oci, OC2) consists of the eigenvectors of 
C22, 

OcC 2 2O' c =( A C V 

Ac < 0, m is the multiplicity of the zero eigenvalue and the columns o/Oci and 
Oc2 consist of eigenvectors with nonzero and zero eigenvalues, respectively. 

Some of the settings described in Lemma 3 concern ratios that are not in the class C-r- 

Lemma 4. When considering the right tail, matrices A and B admit a ratio R in the 
class C-jz only for cases 1, 2(b) or 2(c). When considering both the left and right tails, 
then the class C-jz C\Cc encompasses case 1 and the special setting of case 2(c) in which 
C 22 = 0. 



Lemmas 3 and 4 are most easily understood by means of some simple examples. Con- 
sider an distribution for R. Then, 
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and Ai(r) = — r with Aa(r) = 1. Clearly, this is not in the class C-r nor in Cc- Since 
C22 = 1, a scalar, this is case 2(a). 

Next, consider n — 2 and the least-squares estimate of a lag-1 serial correlation in the 
simplest setting with R ~ e^Ai = £2/^1 ■ Note that this has the Cauchy distribution 
when fi — and the support is (I, r) = (—00,00). To sec that this ratio is in the classes 
Ciz and Cc, note that 

A-r-B=f- r ^ (6) 



1/2 



and that the limiting eigenvalues are 



lim (-r - \Jr 2 + l)/2 = = lim (-r + V ' r 2 + l)/2. (7) 

r — ' — 00 r — >oo 

The example illustrates a case 2(c) ratio in which C12 = 1/2 and C22 = are scalars and 
^V(C22) % iV(Ci2). The same results hold more generally with least-squares estimates of 
serial correlation from regression residuals. 

Lemma 5. Suppose that R has a nondegenerate distribution in the class C-jz, as described 
in Lemma 4, B > and A has rank of at least one. The upper range of support r < 00 
for R, as given in cases 1, 2(b) and 2(c) of Lemma 3, solves A n (r) =0. If r is an interior 
point of the support of R, then the moment generating function of X r is 

M Xr (s)= (n(l-2.sA 4 )- 1 ^ expj^Y^-}, (8) 

convergent in the neighborhood of zero given by 

2A^M <S< 2XW- (9) 
2.2. C.d.f. saddlepoint approximation 

The saddlepoint approximation is based on the cumulant generating function (c.g.f.) for 
X r , given by Kx r {s) — hiMx r {s). The saddlepoint s is the unique root of 

in the range (9). The approximation of Lugannani and Rice (1980) to first order is 

<f>(w) + ^w)^- 1 -it- 1 }, H0^£[X r ], 
Pv 1 (R<r) = l 1 , K£(0) ., n crvl (11) 
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where $(•) and </>(•) denote the distribution and density function of a standard normal 
random variable, respectively, and 

w = sga(s)y/-2K Xr {§), u = s^K Xr (s). (12) 

A second-order c.d.f. approximation is given in Butler and Paolclla (2007). 

2.3. Density saddlepoint approximation 

The saddlepoint density approximation for fii(r), the density of R at r, is derived in 
Butler (2007), Chapter 12.1, or Butler and Paolella (2007) as 

Mr)= , MS) =M Xr (s), (13) 



where s is the same saddlepoint used in the c.d.f. approximation and which solves (10). 
The factor J r (s) is computed from 

J r (s) = tr(I - 2.sA r )- 1 H r + z/ (I - 2sA r )" 1 H r (I - 2sA r )- 1 v r (14) 

with H,,= P r BP^,. A second-order saddlepoint density is given in Butler (2007), page 
383, or Butler and Paolclla (2007). 

Example 6. For matrices A and B in which J?~ Bcta(m/2, (n — m)/2), the saddlepoint 
density in (13) is 

, B(m/2,(n-m)/2) 
B(m/2, (n — m)/2) 

where B is Stirling's approximation for the Beta function B. 

3. Uniformity of the approximations in r 

The relative errors of Lugannani and Rice's approximation in (11) and the density ap- 
proximation in (13) can be shown to be uniform over [0, r) when in the class Cn- These 
results follow as consequences of deriving their finite limiting ratios as r — > r. The limiting 
ratios are derived in Theorems 9, 14 and 15 below. Our approach to computing these lim- 
iting ratios follows that used in Jensen (1988), (1995), Chapter 9.4, and generalizes these 
results to accommodate both nonccntrality and the special concerns involving multiple 
eigenvalues. 

The nature of these asymptotics is dependent on the multiplicity of the eigenvalue 
A n (r) = 0, denoted as m. As a simple eigenvalue with m — 1, the limiting ratios are 
derived in Theorem 9. This is a common setting encountered when dealing with serial 
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correlations. With m > 2, however, the asymptotics are more difficult and such results 
are deferred to Theorem 14. Examples of the multiple eigenvalue setting are also common 
and include least-squares estimates and Yule- Walker estimates for lag-/ serial correlation 
with I > 2. One important multiple eigenvalue example is the noncentral beta distribution 
discussed in Section 4. 

The case 2(a) setting is not in C-r.; however, the relative error can still be shown to be 
uniform over [0, oo); see Butler and Paolella (2007). 

3.1. Simple eigenvalue A n (r) = 

Suppose that r < oo, under the circumstances of cases 1, 2(b) or 2(c). We assume here 
that A — rB has a simple zero eigenvalue with multiplicity m = 1. For general A and B, 
this multiplicity is often difficult to anticipate. Define 



where p„(r) is the eigenvector associated with the zero eigenvalue of A — rB. 
The situation with r = oo is more complicated. 

Lemma 7. Suppose case 2(c) with r = oo. Then, m, the multiplicity of zero eigenvalues 
in {Ai(oo)} ; is the number of zero eigenvalues for C22 ■ If m = 1, then 



where o„ is the px 1 eigenvector associated with the zero eigenvalue of C22, O b = 
(Obi,Ob2) and Ob2 is n x p and the orthonormal basis for the null space of B used to 
determine C 2 2- 

The AR(1) example in (6) with n = 2 provides a simple example. Here, C22 is the 
scalar so that o' n = 1, and B2 = (0, 1). Hence, vq = [J-2- 

Lemma 8. Suppose that the conditions of Lemma 5 hold and let m = 1. Then, as r — > 
r < 00, e = A„(r) — > A n (r) = and s = to/e + 0(1) — > 00, where 



v a = v„(r) :=p„(r)'/i, 



(15) 



v = v n (oo) :=o' t O B2 M, 



(16) 



t = -—{In - l + 



o 2 - y / K 2 + 2n-l) 2 - (2n- 1)2 + 1} 



(17) 



and vq is defined in (15) or (16). In addition, 




(18) 



Theorem 9. Suppose that n>2,R has a nondegenerate distribution in C-r, B > and 
A has rank of at least one. If m = 1, then the limiting ratio of the true tail probability to 
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its first order Lugannani-Rice approximation in (11) is 



PijR>r) _ v /2^(l-2t )(2t )("- 1 )/ 2 Mo e-^ fn 1 
™K 1 ( J R>r)" B(l/2,(n+l)/2)(n/2) ^{2' 2' 2 > 



where 



and parameters io,Mo and are specified in (17), (18) and (15), (16). The first- order 
saddlepoint density has the same relative limit. All of these parameters are determined 
by vq, so the right-hand side of (19) is a function of vq alone. 



In the central case with v = 0, the limiting ratio of tail probabilities in Theorem 9 is 

(i> ^^0/-®(|' wnere B is Stirling's approximation. This same limiting error was 

derived in Jensen (1995), Chapter 9.4, which considered the tail ratio for the distribution 
of the least-squares estimate in a mean zero AR(1) model. Jensen's (9.4.7) is this value 
when the difference in notation is accounted for (our n being n + 1 in that paper) . 

As vl -> oo, the limiting ratio in Theorem 9 is f (|, ^)/r(i, ^){1 + O(^ " 2 )}, where 
r is Stirling's approximation. This follows from the large argument asymptotics for \F\ 
as given in 13.1.4 of Abramowitz and Stegun (1972). 

Relative errors for the second order cdf and density approximations are also uniform 
in the right tail; see Butler and Paolella (2007). 



3.2. Multiple eigenvalue A Tl (r) = 

In this setting, the asymptotics depend on the relative rates of convergence to zero for 
the multiple eigenvalues of A — rB that approach as r — » r. If m denotes its multiplicity, 
then m > 1, by the definition of r. The allowable values of m arc 1 < m < n — 1, but not 
m = n. This latter value would make lim r _» r (A — rB) = 0, in which case the distribution 
of R approaches a degenerate distribution, by Lemma 2. For unbounded ratios in case 
2(c), the value of m is the dimension of the null space for C22, whereas, for ratios in £>, 
the value of m is less transparent. 

We must first determine the relative rates at which the m largest eigenvalues of A — ?'B 
vanish as r — > r in the two separate settings, r < 00 and r = 00. For the former setting, 
general formulae for these relative rates are given in the next lemma. When r = 00, the 
relative rates must be determined on a case-by-case basis. 

Lemma 10. Suppose that r < 00 and is an eigenvalue of multiplicity m for A — rB. Let 
the columns of the nxm matrix Uo be an orthonormal basis for the null space of A — rB . 
Furthermore, denote the ordered eigenvalues of UgBUo as < T„_ m +i < ••• <r n . If 
T n > 0, then the limiting relative rates of convergence to zero for the m largest eigenvalues 
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of A — rB are 



Iim - = — = a>j (21) 
r->rA„(r) r„ 



/or i = ?i — m + 1, . . . , n, where < w„_ m +i < • • • < W„ = 1 . 

For the most common case, in which B > 0, we have T„_ m+ i > so that w n _ m _|_i > 0. 
To deal with the r = oo setting of case 2(c), we reparametrize 

D(e) = (A-rB)/r = eA-B (22) 

and let e = 1/r — > 0. If Aj(er) are the ordered eigenvalues of A — £ _1 B, then 

&(e)=eA<(e) (23) 

are the ordered eigenvalues of (22). 

Lemma 11. Consider case 2(c), in which r = oo, and assume that the zero eigenvalue 
has multiplicity m. A n _ m +i(r), . . . , A„(r) are </ien analytic at r = oo. J/ A^(oo) > 0, then 
the relative rates of convergence are 

_ d\{e)/de\ £=0 d 2 Me)/de% =0 



d\ n (s)/de\e=o d 2 ^ n (e)/de-\ t 



-2 

for i = n — m + l, 



The limiting noncentrality parameters {i>oi : i = n — m + 1, . . . , n} are more difficult to 
determine for m > 2 because they are expressed in terms of the limiting eigenvectors 
associated with the eigenvalues that vanish. In the case r<oo, it is intuitively clear 
and Lancaster (1964) has shown formally that these eigenvectors are smoothly defined as 
1 — ► r. Let P 2r be n x m and consist of the last m columns of P' r which are the eigenvectors 
for the m largest eigenvalues of (A — rB) (which increase in size with column number). 
P27. is then continuous at r = r and the limiting noncentrality parameters are 

.n— m+1 3 * • • j f0,n)' = P 2 rM- ( 25 ) 

In the unbounded setting with r = 00, let the nx m matrix P2 e consist of the eigenvectors 
corresponding to the largest m eigenvalues of D(e) in (22). The limiting noncentrality 
parameters are then given in (25) with P20 = lim e ^oP2e replacing P2 r - In complicated 
practical examples where these computations are not explicit, these limiting eigenvectors 
are best computed numerically by using a small e > 0. 

Example 12. The least-squares estimate of a lag-2 serial correlation with n = 3 and 
zero mean has the form R = 'e\ and leads to the matrix 



D(e)=eA-B= =Q^ 
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where 

'2<Me)/e 2V+(e)/£ 

q;=( o i o ), ^ ± ( £ )=-i±i v / r 

1 1 



The eigenvectors in matrix Q^. have not been normalized as would be needed to use the 
notation P^.. The limits of the eigenvalues are lim e ^ (e),0,ip + (e)} = (—1,0,0) and 
the limiting normed eigenvectors have P' £ — > I3 as e — > 0. Note that di/j+(s) / de\ e= o = 0. 
Also, the eigenvalues of 



'22 



0' B2 AO B2 









are both zero, as discussed in Lemma 11. The limiting rate 

- r A *(r) i, 

W2 - r*£, ATM - -0 d*il> + (e)/de* ~ T/2 " °' 

The limiting noncentrality parameters are 



Expressions for the limiting relative errors depend on {0^,^}. All summations in the 
remainder of this subsection are over S = {n— m + 1, . . . , n}. 

Lemma 13. Suppose that R is in the class C-r and let m be the multiplicity of the zero 
eigenvalue of A — rB. Then, as 1 — ► r, e = A„(r) — ► A„(r) = and s = to/e + O(l) — > 00, 
where to is the unique solution to 







■ F«>i{ + 7 ^-^1 (26) 



m (0, 1/2) im£/i 5 = {n — m + 1, . . . , n}. In addition, 



u—>uq 



Theorem 14. Suppose that n > 2 and the conditions of Lemma 13 hold. Define the 
operator D (X) to be the density of the random variable X evaluated at zero. The limiting 
ratio of the true tail probability for R to its first order Lugannani-Rice approximation in 
(11) is then 
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where the \ 2 terms are independent random variables. Parameters iju and r)2i, for i G S, 
are 



r lH= Ti 7TL 7> V2i 



^{1-2^)' 2(l-2t w i ) 
The limiting ratio for the density approximation is 



ies lies u ° ) 



(29) 



^ f R (r) Wj 

+ E '''■>>'-<>'l>' l l\i h >! l) A VlkX 2 (l, 2T72fe) + VliXl + VljXl ~ — Xn-m \ 

ieSjGS Ikes u ° ). 

where n 3l = (1 - 2t u i )~ 1 , 

= H r = lim r _> r P r BP r and all \ 2 variates in (29) are assumed to be indepen- 
dent. 

In the case m = 1, the results of Theorem 14 reduce to those in Theorem 9. 

4. Noncentral Beta (f , ^p) distribution 

This distribution has 

A _(i m 

A ^ 

and B = I„ so that r —1 and Wj = 1 . This leads to the explicit expression 

to = \ + ^-{{d - m) - y/{9-mf+49n}, 
2 4n 



where 

m 

igS i=l 

Furthermore. 



i n — m „ I to 20 

"j-2 ■ 



«o = i / — =— + 2i? 



2 " °\(l-2to) 2 (l-2t ) 3 
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Theorem 15. For a noncentral Beta(^, ra ~ m ) distribution with min(m,ri — m) > 1, the 
limiting ratio of the true tail probability to its first-order Lugannani-Rice approximation 
is 



where 



yggl - 2t )W2(2t )("-m)/2 Moe -^ /nm L 
B(m/2,(n-m)/2)(n-m)/2 1 ^2' 2 '2/' 



% = 



n — ni \ 



2(1 -2t ) 

XTie first-order saddlepoint density has the same relative error limit. 

In the central setting with = 0, the value in (30) reduces to B(f , ^^)/B(f , 
This is consistent with the computation of the central Beta^, "~ m ) density in Example 
6. 

As 6^oo, the limiting ratio for (30) is f(±, ^)/r(±, ^){1 + O^ 1 )}, which 
follows from the asymptotics for 1F1 given in 13.1.4 of Abramowitz and Stegun (1972). 



5. Examples 

5.1. Serial correlations 

Least-squares, Yule-Walker and Burg estimates for lag-Z correlations are considered in 
further detail in Butler and Paolella (2007) as members of the respective classes C-r—B, 
B and B. Both simple and multiple eigenvalue settings occur in such examples when 
there is no correction for mean. 

In practice, serial correlations with arbitrary lag I are generally computed from least- 
squares residuals and this often ensures that the largest eigenvalue of A — rB has algebraic 
multiplicity one. Thus, the simpler situation for the large deviation errors occurs most 
often in practical data analysis. 



5.2. Numerical example 

Numerical confirmation of the large deviation errors in Theorem 9 is possible by consid- 
ering the simplest model of Section 2.1. This is the least-squares estimate of lag one with 

n = 2 in a model without a location effect. Then, R = 62/61 with ej ln ~ p N(/ij, 1). The 
exact density can be expressed as 

1 f°° , f 1, ,ol f 1 



M r ) = ^ J Mcxpj-^-^) 2 jcxpj--(rx-^ 2 ) 2 \dx 



(n6)- 1 expi--(ni+ f ii)} + 



(31) 



2 V "\ SV2^6 
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-20 -15 -10 -5 5 10 15 20 

Figure 1. Exact density fn (solid), second order fn2 (dashed) and normalized fa (dotted) 
approximations . 



8 = 1 + r\ 9 = erf (tl±^l\ , A = exp( —- ^ 



where 



26 J V 2 5 

From Theorem 9, the limiting relative errors in the left and right tails are dependent on 
vo alone; in both tails, this value is vq = fi2 so that the limiting relative errors are the 
same in both tails, regardless of the values of /Ui and \i2- 

The density (31) is both heavy-tailed and bimodal for [i\ = 0.2 and ^ 2 = 2. Figure 
1 plots the exact density, the normalized version of Jr in (13) denoted by Jr and the 
second-order saddlepoint /_R2 given in (20) of Butler and Paolella (2007). While both 
appear highly accurate in the tails, only the latter captures the bimodality. Figure 2 
plots the ratio of the exact to the three approximate densities including Jr, Jr and Jr2- 
As \r\ increases, we have numerically confirmed that fn{r)/fR(r) — > 0.8222, in agreement 
with the value computed via Theorem 9. This value is virtually achieved at \r\ = 10. Both 
fR and fR2 perform better than Jr in the tails, the latter most notably so. 

The true cdf of R, or Fr(t), must be computed from (31) using numerical integration. 
In this case, |r| must be substantially larger before the same limiting ratio, as specified 
in Theorem 9, is reached. Figure 3 plots 

FrXt) 1 - F R {r) 
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Figure 2. Error ratios /h//h2 (dashed), fa/fa (dashed-dot) and fa/ fa (dotted). 



with Fn(r) as the second-order approximation Pr2 in (17) of Butler and Paolclla (2007) 
and also as Pri in (11). For these values of fii, Pri is more accurate than Pr2 only in 
the range -1.8 < r < 1.2. At r = -25,000, F R (r)/F R (r) = 0.8226 for Pr x , as given by 
Theorem 9, while, for Pr2, the ratio is 1.015. This latter ratio necessarily includes the 
factor (1 + Of): where Of approximates the limit of the second-order correction term. 

If /ii = fi2 = 0, then R is Cauchy and the saddlepoint density reduces to /ij(r) = 
y/-K/2fii{r). Thus, fn is exact and the saddlepoint solution to = K' x (s) is given by 
s = r. The relative error is, B(l/2, l/2)/S(l/2, 1/2), in agreement with the large sample 
space theory. 
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